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Similar to the radiative parton energy loss due to gluon bremsstrahlung, elastic energy loss 
of a parton undergoing multiple scattering in a finite medium is demonstrated to be sensitive to 
interference effect. The interference between amplitudes of elastic scattering via a gluon exchange 
and that of gluon radiation reduces the effective elastic energy loss in a finite medium and gives 
rise to a non-trivial length dependence. The reduction is most significant for a propagation length 
L < 4:/nT in a medium with a temperature T. Though the finite size effect is not significant for 
the average parton propagation in the most central heavy-ion collisions, it will affect the centrality 
dependence of its effect on jet quenching. 
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One of the remarkable phenomena observed in central 
nucleus-nucleus collisions at the Relativistic Heavy-ion 
Collider (RHIC) is jet quenching as manifested in the 
suppression of high transverse momentum hadron spectra 
[1,2], azimuthal angle anisotropy [4] and suppression of 
the away-side two-hadron correlations [5]. The observed 
patterns of jet quenching, their centrality, momentum 
and colliding energy dependence in light hadron spectra 
and correlations are consistent with the picture of par- 
ton energy loss [6]. However, recent experimental data 
on measurements of single non-photonic electron spectra 
[7] seem to indicate a suppression of heavy quarks that 
is not consistent with the theoretical predictions based 
on current implementation of heavy quark energy loss 
[8--10]. This has led to a renewed interest in the elastic 
(or collisional) energy loss and its effects in the observed 
jet quenching phenomena [11-13]. 

Elastic energy loss by a fast propagating parton in gen- 
eral is caused by its elastic scattering with thermal par- 
tons in the medium through one-gluon exchange. The 
gluon exchange can be considered as the emission and 
subsequent absorption of gluons by two scattering par- 
tons and therefore should also be subject to formation 
time of the virtual gluon due to interference. In a finite 
medium when the propagation length is comparable to 
the formation time, the interference should reduce the 
effective elastic energy loss as compared to an infinite 
medium. Since the typical energy exchange with a ther- 
mal parton is w ~ T, the average formation time is there- 
fore controlled by the temperature of the medium. In a 
recent study by Peigne et al. [14], however, collisional en- 
ergy loss was found to be suppressed significantly for con- 
siderably large medium size. This might be partially due 
to the complication of subtraction of induced radiation 
associated with the acceleration of color charges within 
a finite period of time in the semi-classical approach. 

For partons produced via a hard process that go 
through further multiple scattering, elastic scattering via 
one-gluon exchange often bears many similarities to the 



radiative (or gluon bremsstrahlung) processes. In this 
paper, we will consider the elastic scattering within the 
same framework of multiple parton scattering that was 
employed to study the radiative energy loss [15,16]. We 
will show that it is the same interference between the am- 
plitudes of elastic scattering and induced radiation that 
leads to the reduction of the elastic energy loss. Such a 
reduction gives rise to a non-trivial medium size depen- 
dence of the elastic energy loss. 

For the purpose of illustration and derivation, we start 
with the double quark scattering processes in deeply in- 
elastic scattering (DIS) off a nucleus and their effects on 
the nuclear modification of the effective quark fragmen- 
tation functions. The differential cross section of semi- 
inclusive processes e(Li) -I- A{p) — > e(L2) + /^(^/i) + X 
in DIS can be expressed in general as 
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and the leptonic tensor, L^i, ~ ^ Tr(7 ■ Lijf^j ■ L2ji^) , 
where q — [—Q'^/2q^,q^,0±] is the four-momentum of 
the virtual photon, p = [p+,0. Ox] is the momentum of 
the target per nucleon, s = {p + Li)'^ is the total invari- 
ant mass of the lepton- nucleon system, J^ is the hadronic 
electromagnetic (EM) current, J^ = eqipqjf^Tpq, and J2x 
runs over all possible intermediate states. 

The leading-twist contribution to the semi-inclusive 
hadronic tensor to the lowest order in the strong cou- 
pling constant comes from a single virtual photon and 
quark scattering. 
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where ff{x) is the quark distribution in the nucleus, 
X = Q^/2p'^q~ the Bjoken variable, Dq^h{zh) the quark 
fragmentation function and 
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the hard part of 7* + g partonic scattering. 

For the simplest case, let us consider multiple scat- 
tering between two quarks with different flavors in DIS 
off a large nucleus as illustrated in Fig. 1 (with the cen- 
tral cut). In this process, a quark (q) knocked out by 
the virtual photon undergoes a secondary scattering with 
another quark q' from the nucleus. It contributes to the 
semi-inclusive DIS at twist four similarly as multiple scat- 
tering with gluons in the discussion of induced gluon ra- 



diation. We assume the flavors of the two quarks to be 
different so that there is no contribution from crossing 
diagrams. The contribution from the central-cut dia- 
gram in Fig. 1 to the semi- inclusive tensor is, as given 
in Ref. [15], 
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is the two-quark correlation function in the nucleus. The 
four terms with different phase factors in the above 
matrix element correspond to hard-soft, double hard 
quark scattering and their interferences, similarly to the 
Landau-Pomeranchuck-Migdal (LPM) interference in the 
processes of induced gluon radiation [15]. 




FIG. 1. Cut diagrams for double quark-quark scattering in 
DIS off a nucleus. 



In the hard-soft quark scattering, the first or lead- 
ing quark that was knocked out of the nucleus by the 
hard virtual photon becomes off-shell. It then radiates 
a real gluon which interacts with another soft quark 
that carries zero fractional momentum ( + component) 
in the coUinear limit (neglecting the initial transverse 
momentum of the quark) and converts it into the final 
quark with momentum i' . In this case, the final to- 
tal longitudinal fractional momentum of the two quarks 
XL == (^+ + ^'+)/p+ = e^/2p+q-z{l - z) comes from 
the initial leading quark before its interaction with the 
virtual photon. This kinematics corresponds to the in- 
duced emission of a secondary quark like the case of in- 
duced gluon radiation. Note that quarks from the ma- 



trix elements normally are off-shell and carry transverse 
momentum. The kinematics in the hard-soft quark scat- 
tering process will then allow both finite fractional lon- 
gitudinal and transverse momentum when the secondary 
quark is off-shell. It is only in the coUinear approxima- 
tion, i.e., neglecting the initial transverse momentum, 
that the quark lines from the nucleus can be considered 
on-shcll. The kinematics in the hard-soft processes then 
requires the secondary quark to carry zero fractional lon- 
gitudinal momentum. 

In the process of double hard quark scattering, the 
leading quark becomes on-shell after its interaction with 
the virtual photon. It then scatters with another quark 
that carries finite momentum fraction x^. Therefore, the 
final total longitudinal fractional momentum of the two 
quarks xl is transferred from the initial momentum of 
the second quark. This is equivalent to elastic scattering 
between two quarks via one-gluon exchange. 

For a complete calculation of double quark scattering 
processes in DIS, one also has to include interference be- 
tween triple quark scattering and single quark scattering 
as given by the left and right cut diagrams in Fig. 1. 
Their contributions to the semi-inclusive hadronic tensor 
are very similar to the central-cut diagram. 
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with the corresponding four-quark matrix elements 
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Note the negative signs before the phase factors and the 
different time orderings as given by the 9 functions in the 
matrix elements from the interference contributions. In 
addition to the three cut diagrams in Fig. 1, one should 
also include virtual corrections which correspond to dia- 
grams with cut quark lines right after the virtual photon 
interaction. The total contribution to the semi-inclusive 
spectra from double quark scattering should include all 
three cut diagrams and virtual corrections from which 
one can then calculate the nuclear modification to the ef- 
fective fragmentation function due to double quark scat- 
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tering. The virtual corrections make the total contribu- 
tion from double quark scattering to the modified frag- 
mentation functions unitary. However, they do not con- 
tribute to the quark energy loss we discuss in this paper. 
To compute the energy loss experienced by the leading 
quark due to double quark scattering, one can identify it 
with the energy fraction the leading quark loses to the 
second quark in the central-cut diagram or the radiated 
gluon in the left and right cut diagrams. When summing 
these contributions, one encounters the combination of 
the three 6*- functions, 
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which is a path-ordered integral limited by the value of 
77^ '^ 1/xp^. We will neglect any such contact con- 
tributions as compared to the dominant contributions 
whose integration over the location of the second quark 
is only limited by the size of the nucleus. In this case the 
two quarks come from two different nucleons in the nu- 
cleus and the corresponding contributions are enhanced 
by A^/^ due to the large nuclear size as compared to the 



contact contribution. The dominant quark energy loss 
due to double quark scattering is then 
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where the four-quark matrix element Tf^,{x,X]^) is de- 
fined as 
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One can notice that the above matrix element has an 
overall phase factor e"^P '^^2 ^Vi ) as a consequence of 
the cancellation by the interference processes in the left 
and right cut diagrams. Therefore, the remaining contri- 
bution from the double quark scattering requires the sec- 
ond quark to carry initial longitudinal momentum xlP'^ 
that provides the total longitudinal momentum of the 
final quarks after the elastic scattering between the lead- 
ing and the second quark. The above contributions also 
contain the interference between elastic and radiative am- 
plitudes. The leading quark in the radiative process first 
emits a real gluon which converts a soft quark (which is 
close to real for small initial transverse momentum ) into 
the final quark with momentum £' . 

To a good approximation, especially in the case of par- 
ton propagation in dense medium in heavy- ion collisions, 
we can neglect the correlation between the initial leading 
quark and the secondary quark in the nuclear medium. 



The above matrix element can be factorized as a product 
of the leading quark f^ (x) and the secondary quark dis- 



tribution in the medium, 
medium, we have 



In the rest frame of the nuclear 
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where y — [yi + y2)/2, fq'{xL) is the secondary quark 
distribution per nucleon and p{y) is the spatial profile of 
the nucleus. 

One can extend the above calculation to quark prop- 
agation in a hot medium in high-energy heavy-ion colli- 
sions. For scattering between a leading quark with energy 
E and a thermal quark with energy u>, the quark distri- 
bution in a thermal medium with quarks as the basic 
constituents at temperature T is, 
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where xl = i'^/2u!Ez{l — z), dq — AN^Nf is the num- 
ber of degrees of degeneracy for the quark and the spa- 
tial dependence is exphcitly through the temperature 
T{y). One can complete the integration over z and ^t in 
Eq. (11). The corresponding elastic energy loss is then 
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where x^^ = -^/l — Ip? jioE and \i is the transverse mo- 
mentum cut-ofF which regularizes the collinear divergence 
normally associated with parton radiation or small angle 
scattering. In vacuum, such a cut-ofF separates perturba- 
tive and non-perturbative regime of QCD. In a medium 
at finite temperature, a resummation of hard thermal 
loops in the gluon propagator can also provide regular- 
ization of the collinear divergence as in previous studies 
[17]. 

With the above equation, one can calculate the elas- 
tic energy loss for any spatial profile r(y) of the thermal 
medium. To simplify the calculation, we neglect the ther- 
mal fluctuation of the quark energy in the variable Xu^ 
and replace lo with ST. One can carry out the thermal 
averaging and obtains, 
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and Pq{T) — dqri{3)T^ /tt'^ is the local thermal quark den- 
sity and XT = y/T~2i?j3TE. 

One can similarly consider quark energy loss due 
to elastic quark-gluon scattering which is more in- 
volved. In particular, the quark-gluon Compton scat- 
tering processes are very similar to the induced gluon 
bremsstrahlung as considered in previous studies of ra- 
diative energy. As we will show in a separate study [18], 
the original result [15] for radiative parton energy loss ac- 
tually contains both radiative and elastic processes. The 
LPM interference is in fact the interference between the 
radiative and elastic scattering amplitudes that are also 
responsible for the interference effects in the elastic en- 
ergy loss in this paper. The final result for the elastic 
energy loss due to quark-gluon scattering is 
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where 
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Pg{T) — dg({3)T^ /tt'^ is the local thermal gluon density 
and dg ~ 2{N^ — 1) is the number of degrees of degeneracy 
for gluons. Note that 77(2) = C(2)/2, r]{3) = (3/4)C(3), 
C(2) = TrVe and C(3) « 1.2021. 

For the simplest parton density profile, we consider 
a finite medium with a constant parton density and a 
length L. The spatial integrations in the above elastic 
energy loss can be carried out. The total elastic energy 
losses are {dE/dL = AE/L) 
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For a medium with infinite length L ^ 00, the above 
elastic energy losses have finite asymptotic values, 
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which can be interpreted as the average energy transfer 
V = l3pl2ijj per mean-free path \qa = ^j Pa^^qa [19] due to 
elastic scattering between the leading quark and thermal 
partons. Here daqajdGp is the differential cross section 
for g + a elastic scattering. 

According to Eqs. (20) and (21), 



in the hmit SET j \? > 1. The extra factor 3/2 in 
Eq. (24) as compared to early result [17] comes from the 
full splitting function [1 + (1 — z')^\l z we used in the cal- 
culation. If we take the approximation z <^ \ which 
corresponds to small angle quark-gluon scattering in the 
early calculation [17], the factor 3/2 will disappear. 



One can observe from the above results that the ef- 
fective elastic energy loss for a quark propagating in a 
medium with a finite length L is reduced from its asymp- 
totic value due to the interference between the elastic and 
radiative amplitudes. The underlying physics is very sim- 
ilar to the LPM interference in the gluon bremsstrahlung 
induced by multiple scattering. In the elastic scattering, 
the gluon exchange can be considered as the gluonic field 
built up around the propagating quark over finite period 
of time. The formation time in the case of elastic scat- 
tering with thermal partons is controlled by the thermal 
energy of the parton and therefore is determined by the 
temperature on the average. Shown in Fig. 2 is the quark 
elastic energy loss dEqa/dL{a — q, g) due to quark-quark 
and quark-gluon scattering normalized to the asymptotic 
value (dEqa/dL) 00(0, = q,g) in a finite medium as a func- 
tion of a dimensionless variable -kTL. The dependence on 
L is the strongest when -kTL <^ 1, 
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When i ^ 0, the destructive interference becomes com- 
plete and the elastic energy loss vanishes. 



expanding system with a realistic parton density profile, 
one has to numerically compute the elastic energy loss 
and the corresponding modified fragmentation functions 
according to Eq. (15). Furthermore, one should also treat 
radiative and clastic energy loss in the same framework 
since there is intricate connection between the two. This 
will be discussed in a separate work [18]. 

Note added: After completion of this paper, the au- 
thor noticed a recent publication [20] in which the finite 
size effect on clastic energy loss is a result of the imposed 
finite interaction time which modifies the energy conser- 
vation and therefore the phase space available for final 
partons after elastic scattering. The interference effect 
discussed in the current paper is not included and there- 
fore the obtained length and energy dependence of the 
elastic energy loss are different. 
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FIG. 2. Elastic energy loss in a finite size thermal medium 

(at temperature T) dE/dL as a function of the medium length 

L normalized to the asymptotic value {dE/dL)oc for L = 00. 

The typical average propagation length of a parton 
in the most central Au + Au collisions is about L ^ 6 
fm. For an average temperature T ^ 200 — 300 MeV, 
TiTL ^ 19 — 28. In this case the interference effect is quite 
small as seen in Fig. 2. This is quite different from the re- 
sults in Ref. [14]. However, for partons produced in the 
peripheral region of heavy-ion collisions whose propaga- 
tion length is on the order of 1-2 fm, the reduction of the 
elastic energy loss due to interference is still significant. 
For an accurate and consistent treatment of elastic en- 
ergy loss, one should take into account the interference 
effect and the non-trivial distance dependence. For an 
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